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Abstract 



We have analysed the recent results from the observation of charged cur- 
rent v e d — > e~pp events from solar neutrinos by the Sudbury Neutrino Ob- 
servatory SNO assuming neutrino oscillations with three active flavours. The 
data seem to prefer a low mass-squared difference and large mixing angle so- 
lution (the so-called LOW solution) in (12) parameter space. However, when 
combined with the Gallium charged current interaction data from Gallex and 
GNO, distinct (la) allowed regions corresponding to the large mixing angle 
(LMA) and small mixing angle (SMA) appear while the LOW solution is dis- 
favoured upto 3a standard deviation. The physical electron neutrino survival 
probability corresponding to these best fit solutions are then determined and 
analysed for their energy dependence. 
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I. INTRODUCTION 



The recent announcement of results on the observation of charged current (CC) neutrino 
events from the sun by the Sudbury Neutrino Observatory (SNO) |1] has lead to a flurry of 
activity in understanding the implications of the data on neutrino oscillations. Several au- 
thors have concentrated on the global analysis of the solar neutrino data from all existing 
observations and the refinement in the allowed oscillation parameter space. Typically, such 
analyses have refined the already existing allowed regions in parameter space — the so-called 
small mixing angle (SMA), large mixing angle (LMA), Low Mass (LOW) and the 'just so' 
solutions. However, not as much attention has been paid to how individual experimental 
results perform with respect to the best possible global fits. 

This paper is not devoted to providing yet another global best fit. We concentrate rather 
on the physical quantity of interest, viz., the neutrino survival probability. Since we restrict 
ourselves to the CC data, we focus on the electron neutrino survival probability, P ee (E u ). 
The question that we seek to address is, what are the constraints from data on the energy 
dependence of P ee . In order to answer this question, we use the energy integrated CC data 
from the combined Gallium experiments of GALLEX and GNO[] Q (henceforth referred to 
as the Ga experiments) as well as the recent CC observations from SNO |]]. We also briefly 
discuss the constraints arising from the Homestake Chlorine experiment || (henceforth called 
CI). We do not include the results from Super-Kamiokande in our analysis since the observed 
elastic scattering process has both CC and neutral current (NC) components. We use the 
standard solar model of BP2001 in our theoretical analysis; the flux data are available in 
tables at the web-site 0. 

Of these experiments, the Ga experiments observe neutrinos over the widest range of 
energies. However, the Ga data contain no energy information. Theoretically, the bulk of 
the events may be attributed to electron neutrinos with energy less than 1 MeV with less than 
10% coming from neutrinos with energy more than 5 MeV. On the other hand SNO observes 
only the 8 B neutrinos above a threshold of about 8 MeV. It is also sensitive (indirectly) to 
the neutrino energy. In short, these experiments are sensitive to different source reactions 
in the sun. We can therefore get good information by analysing the combined data as well 
as the individual data sets. The CI experiment is sensitive to the (second) 7 Be line apart 
from the B neutrino flux. Like Ga, it does not contain any energy information. However, it 
turns out that CI and SNO yield similar likelihood curves in the allowed parameter space. 

We first analyse the experiments individually and then the combined Ga and SNO data. 
We have used simple expressions for the likelihood contours where the x 2 is defined using 
diagonal errors only; all correlated errors are ignored. This approximation is obviously valid 
for SNO since there is exactly one contributing source flux; it remains true to good accuracy 
for Ga since the dominant low energy pp source flux has very small errors. Finally, since 
there is only 10% overlap between the Ga and SNO source fluxes, the neglect of correlated 
errors is valid while analysing the combined data as well. This is not true for the combined 
Ga, CI and SNO data sets; hence we have not presented the results of such a combined 



x We have used the cited average of these data; the SAGE Gallium experiment H is also consistent 
with these data. 
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analysis here. 

Assuming that neutrino oscillations are the source of the discrepancy between data and 
theory, we find that the likelihood contours in the allowed parameter space for the Ga ex- 
periments differ substantially from the SNO and CI ones. Consequently, it turns out that 
the electron neutrino survival probability (the CC analysis depends only on this probability) 
corresponding to the best fit parameters (x 2 minimum) has a very different energy depen- 
dence for the fits coming from analysing the Ga or SNO data alone and the combined fits. 
This is true especially for the so-called LMA solution where the SNO prefers a rather flat 
P ee while Ga shows a distinct energy dependence. This leads to the question whether the 
low energy pp neutrinos behave very differently from the higher energy ones. If so, it is 
possible that a measurement of the (second) 7 Be line spectrum will be able to settle this 
issue. For example, an almost vanishing 7 Be flux (< 10-20% of the expected flux) will point 
unambiguously to the SMA as the preferred solution to the neutrino oscillation problem. A 
similar analysis for the best-fit global solutions may be found in Ref. [0]. 

In Section 2, we present the relevant formulae for the event rate as well as for the 
survival probability of (e-flavour) neutrinos in the 3-flavour oscillation scheme. In Section 3 
we present the results of the numerical calculation using the standard solar model of BP2001 
H for the solar neutrino fluxes. Section 4 contains the summary and discussion. 

II. RELEVANT FORMULAE 

We briefly list the relevant formulae in a 3-flavour oscillation scenario to explain the 
discrepancy between data and theory. The mixing is expressed as, 

( v e v T ) T = V x ( v\ v 2 zaj ) T , 

where the mass eigenstates V{ have masses m; and the mixing matrix is parametrised in 
terms of three angles^ as a set of 3 independent 2x2 rotations: 

V = V 23 (iP) x Vutf) x V 12 (u) . 

The angle ip does not occur in expressions for solar neutrinos since only u e neutrinos are 
produced in the sun. 

Since we are interested only in the CC events, the only relevant probability is the electron 
neutrino survival probability, P ee . This can be expressed in terms of the (12) and (13) mixing 
angles, u> and 0, as well as the (12) mass squared difference, m 2 . — m\ = <5 12 . We have 

P ee = sin 2 (f) sin 2 (j) m + cos 2 cj> cos 2 <j>. 

Here the subscript m includes the usual (solar) matter effects. We do not include earth 
matter effects here; this affects only the 'LOW solution and that too, at a small level. 
Here Plz is the non-adiabatic jump probability f| and is relevant near resonance. While 



2 We ignore the CP violating phase here. 



- + (- - P LZ ) cos 2^ cos 2o.'„ 



(1) 
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cos20 m ~ cos 20, cos2u; m depends on the vacuum mixing angles as well as 5u and the 
matter term A(E U ). Hence it is energy dependent. Note that we have used the phase 
averaged expression for P ee ; this is allowed as long as we restrict ourselves to 5 i2 > 1CT 8 
eV 2 . Hence we do not consider the 'just-so' solutions in this paper. Finally, A(E U ) is also 
density dependent ||. The solar density varies with distance; however, we have used average 
values of density corresponding to the peak production of different source fluxes since we 
have averaged over the production region. The results are not sensitive to this approximation 
i- 

The survival probability appears in the expression for the ratio: 



( dE v P ee {E v ) <f>(E v ) a{E v ) 

R th = 1 r , (2) 

dE„<f>(E„)<T{E v ) 



where the factor (P ee <ft) in the numerator represents the depleted flux due to mixing and 
depends on the mixing parameters as shown; a is the cross-section for the relevant process 
in the detector; a may also be differential with respect to one or more variables as in the 
case of the SNO data where we use da/dE e , with E e the scattered electron energy. There 
is a further smearing out of this data due to the resolution function since the measured 
electron energy E e is a gaussian distributed around the true electron energy E' e ]I| . The Ga 
cross-section is taken from tables available in Ref |J. The deuterium cross-section is taken 
from the tables available in Ref. Jiy . 



The theoretical value of R th calculated this way is then compared with the ratio R cxp 
measured by the different experiments^. To avoid double counting, the theoretical errors 
were factored into R th while the experimental numbers included only the statistical errors. 

The significance of SNO and Gallex for mass and mixing parameters is obtained by 
separately determining the minimum x 2 , where we use the standard definition 

X 2 = E — 2 j , (3) 

a "a 

where the explicit parameter dependence of the ratio is i?„ = R t ^-(Sij,6ij) for a given set 
of mass squared differences, 6ij, and mixing angles, %. The range of the subscripts is 
determined by the number of neutrino flavours. For P ee we have the three parameters 5u, 
u and cf). The subscript a runs over the set of data points, R^ p , of a given experiment (s). 
The theoretical and experimental errors are added in quadrature in a 2 . Likelihood contours 
are drawn using 

£ = exp[- X 2 /2 • 

While £ by itself has no meaning (it is the log likelihood or y 2 itself that is relevant), na 
likelihood contours can be drawn corresponding to \ 2 — Xmin + ^ 2 - These contours have the 



3 The SNO ratio is calculated using the same BP2001 model. The experimental event rate in 
SNUs for the Ga (actually the average from Gallex and GNO) and CI experiments were used and 
converted to a ratio using the BP2001 calculated prediction. 
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usual probability interpretation of allowed regions. We use the parameters corresponding to 
the best fit (and within la of the best fit) to find the preferred survival probability. 

As we have observed, all charged current events are directly proportional to the electron 
neutrino survival probability. This is not true for the case of the Super-Kamiokande data in 
the presence of an additional sterile flavour; analysis of this data would require an additional 
unknown, the oscillation of u e into the sterile flavour. This is why an analysis of charged 
current data provides the most restrictive bounds on the mixing. 



III. NUMERICAL ANALYSIS AND RESULTS 
A. Choice of data 

We have individually analysed the following sets of solar neutrino data: CC data from 
SNO (referred to as SNO data), from the Gallium detectors GALLEX and GNO (referred 
to as Ga data), and from the Chlorine detector at Homestake (referred to as Chlorine data). 
All these experiments measure charged current data (we have not analysed the SNO elastic 
scattering data here) and hence are particularly tuned to the electron-type neutrinos. The 
data are expressed as ratios of observed to expected event rates. We use the BP2001 
standard solar model for the theoretical calculation. 

Energy integrated data from Ga and CI (with E u ^ min = 0.24 and 0.814 MeV respectively) 
yield one data point from each experiment while SNO (apart from an energy integrated 
value for the ratio of observed to expected rate) also gives data in 11 bins for the recoil 
electron energy starting from E e of about 7.25 MeV, with bin size about 0.51 MeV (except 
for the last bin) The Q value for the reaction v e d — > epp is 0.93 MeV, so the parent 
neutrino energy is at least 8 MeV; it is however useful to note that the cross-section for this 



reaction peaks typically at E e — E v — 1.4 MeV [K| although the peak value is mildly energy 
dependent. Hence there is good energy information in the SNO data. 

As pointed out earlier, SNO sees only the 8 B neutrinos from the sun, Chlorine also sees 
most (about 90%) of the Be neutrinos and Ga sees both Be and B and much of the pp 
neutrinos. In our analysis, we have ignored errors due to correlations between the different 
fluxes (that would not affect SNO) contributing to the event rate at any given experiment. 
We have also used only the statistical errors in the data (since, typically, systematic errors 
are also correlated between data bins). Hence the x 2 minima that we compute are a lower 
bound on that which would have been obtained if the correlated errors were taken into 
account. 

Also, as pointed out in the Introduction, it is possible to do a combined analysis of the 
Ga and SNO data while still ignoring correlated errors. This has also been done. 



B. Allowed Parameter Space 

We begin by finding the allowed parameter space for the individual data sets. These are 
shown in Figs. [1] where the allowed la regions of 5\2 and u around the global minimum \ 2 
are shown for = 0. There is not much sensitivity to in the range allowed by CHOOZ, 



< 9° [ITJ. The la regions (not shown in the figure) are quite large, especially for the Ga 
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and CI data sets since there is exactly one data point (total rate) available. While it may 
appear that one data point is being fitted with two (or even three) free parameters for Ga 
and CI, the resulting restriction in parameter space occurs because of the non-trivial energy 
dependence of P ee and the fact that different source fluxes contribute differently at different 
energies. For the SNO data, we have shown the region of parameter space allowed by taking 
into account either all the 11 electron energy bins or the integrated rate. Obviously, the 
allowed region is larger if only the average (integrated) rate is taken into account, as shown 
in the figure. However there is a good overlap between these two fits, indicating the stability 
of the fits. 

There is almost complete overlap between the allowed parameter space from Chlorine 
and SNO data especially at large angles. The contribution to CI from the region below the 
SNO threshold comes roughly equally from the (second) 7 Be line spectrum and from a part 
of the 8 B flux. The overlap between the two sets therefore seems to indicate that the nature 
of suppression in the SNO energy regime is similar to that at lower energies, down to the 
CI threshold. The SNO energy-dependent data are more restrictive; the effect of this is to 
allow a narrower strip of parameters than the Chlorine data. 

What is interesting is that the preferred regions for Ga and Chlorine/SNO are totally 
different as can be seen by comparing the relevant panels of Fig. [TJ. The so-called LMA 
solution for Ga has a well-defined u ~ 25° but accomodates a large range of <5i 2 , £12 > 10~ 5 
eV 2 . Solutions corresponding to the SMA exist although the SMA region is not distinct 
from the LMA one; however, the LOW solution (corresponding to low mass squared, large 
angle) is conspicuously absent in Ga. 

Even more surprisingly, while SNO shows allowed regions in the regions corresponding 
to SMA, LMA and LOW, these regions are connected with no well-defined LMA or SMA 
minima. The global minimum for the entire region is in fact the so-called LOW solution 
centred around 5u ~ 6.7 x 10 -8 eV 2 and 00 ~ 29°. 

In general, the SMA and LMA allowed regions of SNO and Ga are different. Hence, when 
the Ga and SNO data are combined distinct regions in LMA and SMA parts of the parameter 
space appear (see Table 1). Since Ga does not have a LOW solution, the 'LOW solution 
(that is in fact the global minimum of the fits to the SNO data) is completely disallowed 
by the combination of Ga and SNO data by about 4cr. The SMA solution is marginally 
preferred over the LMA. However, in our simple analysis we have neglected correlations as 
well as earth matter effects. Though small, this precludes us from preferring SMA over 
LMA; within the errors of the analysis both are equally likely. 

The result from the combined fit is shown in Fig. |[ Note that since we have ignored 
correlations between the fluxes, the allowed regions are just the overlaps of the separately 
allowed regions. In the figure, the la allowed contours (with respect to the local \ 2 minima) 
are shown in the 5u-uj plane for fixed <fi = 0°, 9°, and 20°. The solid lines correspond to the 
combined analysis of Ga and SNO energy-binned data, while the dotted lines correspond to 
the combined Ga and integrated SNO data. As expected, the allowed region shrinks when 
the energy dependence of the SNO data is taken into account. We have also analysed the 
Ga data together with data in the higher energy bins of SNO (from the 5th to the 11th bin 
in the recoil electron energy). This is because the correlated errors are significant (compared 
to the statistical errors) amongst the first four bins and negligible for the rest. Hence our 
our neglect of correlated errors is strictly valid only for these bins. The result of such an 
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analysis (for = 0) is very similar to the analysis using the entire SNO data set and is 
shown in Fig. ^| as dashed lines in the second panel. This leads us to believe that our results 
are not very sensitive to inclusion of correlations. 

Also, for the LMA solution, the \ 2 is smallest for = 0°, increasing with increasing 
4>; the opposite behaviour is seen for the SMA solution. In both cases, however, there is a 
decrease in the allowed region with increasing cf). This provides a loose bound on <ft, <ft < 20°, 
independently of chooz data. This is also seen from the global minimisation with respect 
to all three free parameters. The LMA best fit yields <p ~ with a la error of about 23° as 
can be seen from Table 1. 

C. The survival probability 

The probability shapes corresponding to these preferred solutions are shown in Figs. |3| 
and [J[ The two figures correspond to the so-called large mixing angle (LMA and LOW) 
solutions and the small mixing angle (SMA) solutions. While all three solutions are allowed 
only by SNO, Ga has a preferred SMA and LMA solution, whose corresponding P ee is 
shown in the figures. The P ee corresponding to the \ 2 minimum from the combined Ga- 
SNO analysis is shown as solid lines in the figures. The two lines per solution correspond to 
the la limits on the parameter values and indicate the range of P ee allowed. 

The energy corresponding to the (second) 7 Be line is indicated on both the figures for 
LMA/LOW and SMA solutions. The LMA and SMA solutions clearly show that roughly 
60% and 10% respectively of the total Be flux will survive and be detected at Borexino, 
for example. Even beyond the la level, a low (less than about 30%) fraction of 7 Be will 
indicate a clear preference for the SMA solution. Hence a measurement of the Be flux will 
be crucial in distinguishing the LMA and SMA solutions and in narrowing down the allowed 
parameter space for mixing. 

It is seen that the results of the combined analysis are driven mostly by the SNO data. 
Hence the errors on the P ee are smaller at larger energies where SNO data are available. 
In particular, while still being well within the range allowed by the Ga SMA solution, the 
allowed P ee range from the SNO SMA solution at larger energies is very restricted as can be 
seen from Fig. f|. 

The fits favoured by Ga data, on the other hand, tend to either underestimate (Ga-LMA) 
or allow a very broad band of P ee in the energy range relevant for SNO. The reason for this 
is obvious: there is no energy information in the Ga data. In spite of this, it is able to rule 
out the LOW solution for SNO. (On the average, the Ga data corresponds to lower energy 
of roughly (E u ) ~ 0.3 MeV. The observed Ga ratio which is 0.58 ± 0.07 cannot therefore be 
accomodated by the LOW solution as can be seen from Fig. |3[). 

Another way of understanding this is to look at the likelihood contours for the combined 
Ga-SNO fits over the entire parameter region. We have shown this in Fig. [| for <fi = 0. The 
likelihood contours corresponding to na allowed regions, n = 1, 10, with respect to the global 
minimum (see Table 1) are shown in the 5\2 vs tan 2 a; (or u) plane. The la contours show 
the two distinct SMA and LMA allowed regions very clearly. It is seen that the minimum 
coresponding to SMA is very narrow. The LMA region is larger and smoothly goes into the 
so-called "LOW" at 4a. Hence the LOW solution is allowed only at 4a. The results are not 
very different for = 9°. 
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IV. SUMMARY 



We have performed an analysis of the charged current SNO and Gallium solar neutrino 
data with a view to understanding the physically allowed electron neutrino survival proba- 
bility assuming an oscillation hypothesis with 3 active neutrino flavours. We have analysed 
the data separately as well as done a combined fit. The results can be summarised as follows. 

• The data are not very sensitive to the (13) mixing angle, <fi. However even if we do 
not apply the CHOOZ bound [O], < 9°, very large > 20° does not seem compatible 



with the SNO/Ga data. 

The individual Gallium and SNO (as also the Chlorine data) give rise to very different 
allowed regions of parameter space. These regions are contiguous with no distinct 
LMA and SMA regions visible (see Fig. [T]. 

In particular the so-called 'LOW solution (small (12) mass-squared difference, 5 12, 
large (12) mixing angle, uj) is totally incompatible with the Ga data although it is the 
preferred solution for SNO. A combined analysis of the SNO and Ga data rules out 
this solution at 4(7. 

Distinct small angle and large angle solutions (SMA and LMA) appear on combining 
the Ga and SNO data sets. Both solutions are equally preferred. 

Finally, the currently allowed SMA and LMA solutions predict very different survival 
probability of the (second) 7 Be line spectrum. An exclusive measurement of the 7 Be 
flux such as from Borexino will be able to distinguish whether the small angle solution 
or the large angle solution is ultimately the right solution. 



Due to the severity of the CHOOZ constraint fl2|l , it is possible that these results are not 
significantly altered in the presence of a fourth (sterile) neutrino flavour. 

It may be pointed out that other analyses find all three allowed regions, LMA, SMA 
and LOW. While the Ga experiments yield one data point (since it is an energy integrated 
measurement), the SNO and Super-Kamiokande experiments correspond to many points 
because of energy information. It is possible that the significance of the Ga data is reduced 
when doing a combined y 2 analysis of such a data set. 

A combination of the Ga, CI and SNO (integrated) data can be used to determine certain 
"energy-averaged' values of P ee [0]. For example, a constant P ee corresponding to the SNO 
data at an average energy (E v ) = 10.8 MeV may be assumed. Then the SNO data determine 
the value of P ee at this energy. Assuming that this value holds for all 8 B neutrinos, we obtain 
the central value, P ee (E®) = 0.347 from SNO data. The CI data has contributions from the 
8 B flux as well as the 7 Be flux. Using P ee (E®), we obtain P ee (E®°) = 0.29; finally, using 
these values in the Ga data, the average survival probability in the pp region is obtained to 
be P ee (EP p ) = 0.76. Since (E™) = 0.3 MeV, and E* e = 0.86 MeV, while (E*) ~ 9 MeV, 
it therefore appears as though there is a large energy dependence of P ee in the lower energy 
E v < 1 MeV regime. However, it is crucial to realise that this simplistic analysis may not 
hold when a detailed energy analysis is considered. For example the LMA solution that 
we have obtained corresponds to a flatter P ee around 0.6 in this low energy region while 
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the SMA corresponds to a steeper one with an almost vanishing P ee (E®°. In either case the 
energy averaged predictions do not hold in detail. A completely unambiguous determination 
of P ee in the E v < 1 MeV regime must therefore await either results from Borexino or some 
new energy-sensitive data in the pp sector. 
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FIGURES 




FIG. 1. la allowed contours at <j) = for the individual Chlorine, Gallium and SNO data sets. 
Results for the SNO total integrated data as well as for the energy-binned data are separately 
shown in the bottom panels. 
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FIG. 2. Allowed la regions in 612 an d the (12) mixing angle uj (both uj and tan 2 to are shown) 
parameter space for the combined Gallium and SNO data sets. The solid lines use the energy-binned 
SNO data while the dotted lines use the integrated SNO data. Allowed regions are shown for fixed 
values of 4> = 0, 9, 20 degrees. It is seen that the allowed region shrinks with increasing <p. The top 
right panel shows the fits (at = 0°) obtained by using all the SNO binned data (solid lines) and 
that from excluding the first four bins where correlated errors are large (dashed lines). 
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FIG. 3. The electron neutrino survival probability P ee as a function of neutrino energy for the 
best fit solution in the LMA/LOW regions for = 0°. Two sets of lines for each solution indicate 
the la band in the parameters. Solutions to individual SNO and Ga data are shown as dotted and 
long-dashed lines; the SNO LOW solution is shown as dashed lines. The solid lines indicate the 
preferred P ee from fits to the combined Gallium and SNO data sets. 
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FIG. 4. The electron neutrino survival probability P ee as a function of neutrino energy for the 
best fit solution in the SMA region for <f> = 0°. Two sets of lines for each solution indicate the 
la band in the parameters. Solutions to individual SNO and Ga data are shown as dotted and 
long-dashed lines. The solid lines indicate the preferred P ee from fits to the combined Gallium and 
SNO data sets. 
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FIG. 5. Likelihood contours in the (5i2-tan 2 w plane for <ft = 0° corresponding to na deviations, 
n = 1, 10, from the global minimum value of x 2 fo r the combined Ga and SNO data fits. The LMA 
and SMA allowed regions are marked. The central regions of the parameter space are excluded by 
more than 10a. 
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TABLES 



A, „ 
012 




Ijj 


9 


X 


(8.5(-5.9,+2.8)) 


x 10- 5 


32.0(-4.9,+4.7) 





3.1 


(5.9(-3.1,+2.1)) 


x 10- 6 


2.5(-0.4,+1.2) 





2.6 


(1.0(-0.4,+0.2)) 


x 10- 4 


30.2(-4.8,+5.6) 


9 


3.3 


(7.2(-1.6,+0.7)) 


x 


2.2(-0.2,+1.3) 


9 


2.6 


(1.2(-0.1,+0.3)) 


x 10" 4 


23.9(-4.8,+5.0) 


20 


3.9 


(7.3(-0.4,+0.5)) 


x i(r 6 


2.0(-0.1,+0.1) 


20 


2.4 


(8.7(-6.1,+4.8)) 


x 10" 5 


32.1(-10.9,+4.9) 


0.003(-23.0,+23.0) 


2.8 



TABLE I. Values of 5u and oj corresponding to x 2 minima for the combined fits to SNO and 
Ga data for different fixed values of (f>. la errors on the parameter values are also given. The last 
row corresponds to a global minimum with respect to all three parameters, S±2, 10, and (p. 
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